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Abstract 

We present a derivation of the stress field for an interacting quantum system 
within the framework of local density functional theory. The formulation is 
geometric in nature and exploits the relationship between the strain tensor 
field and Riemannian metric tensor field. The resultant expression obtained 
for the stress field is gauge-invariant with respect to choice of energy den- 
sity, and therefore provides a unique, well-defined quantity. To illustrate this 
formalism, we compute the pressure field for two phases of solid molecular 
hydrogen. 
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The stress, or the energetic response to deformation or strain, plays an important role 
in linking the physical properties of a material (e.g. strength, toughness) with the behavior 
of its microstructure. In addition, the spatial distribution of stress is an invaluable tool 
for continuum modeling of the response of materials. The stress concept has been applied 
at atomistic scales as well. Over the last fifteen years, there has been a continuing trend 
toward understanding various structural and quantum-mechanical phenomena in materials 
in terms of their response to stress 

For example, the residual stress at equilibrium has been used to assess the structural 
stability of systems containing surfaces or strained interfaces. It has been demonstrated 
that the desire to minimize surface stress can give rise to reconstructions on high symmetry 
surfaces and the stability of epitaxially grown bimetallic systems has been attributed 

to the formation of incommensurate overlayers, defects, and dislocations which minimize 
the stress near the metal- metal interface P,|10[. The stress can have significant effects on 
chemical reactivity as well. It has been shown that small molecule chemisorption energies 
and reaction barriers on certain strained metal and strained semiconductor surfaces are quite 
different from those on the unstrained surface 1TT|.[T^ . 



Formally, studies of the above phenomena must include a quantum-mechanical descrip- 
tion of the system's electronic degrees of freedom. Therefore, one must consider how a stress 
is defined quantum mechanically. Methods for calculating the stress in quantum mechan- 
ical systems have been developed since the birth of quantum theory itself [[HJ. However, 
research in developing formalisms for determining the quantum stress in solid-state systems 
has recently been revitalized. This is mainly due to ever-increasing opportunities to per- 
form accurate and efficient quantum-mechanical calculations on systems which exhibit stress 
mediated phenomena. 

The stress is a rank-two tensor quantity, usually taken to be symmetric and therefore 
torque-free. Two useful representations of the stress tensor are the volume-averaged or total 
stress, T a/3 , and the spatially varying stress field cr a( g(x). The two representations are related 
since the total stress for a particular region in a system is the stress field integrated over 



the volume. Nielsen and Martin developed a formalism for calculating the total quantum 



stress in periodic systems [Il|. They define the total stress as the variation of the total 
ground-state energy with respect to a uniform scaling of the entire system. This uniform 
scaling corresponds to a homogeneous or averaged strain over the entire system. They 
further demonstrate that the total quantum stress is a unique and well-defined physical 
quantity. Their formulation has been successfully implemented to study a variety of solid 
state systems [@|§. Other formalisms for determining the total quantum stress have been 



created as well 12.15. 16 



Although these formalisms have provided important tools for studying quantum stress, 
the stress field is a more useful quantity that contains important information regarding the 
distribution of the stress throughout the system. A knowledge of the spatial dependence of 
the quantum stress is vital if one wishes to predict the spatial extent of structural modifica- 
tions or understand phenomena at interfaces in complex heterogeneous systems. However, 
certain definitions of the quantum stress field suggest that it can only be specified up to a 
gauge. (This ambiguity manifests itself in classical atomistic models as well.) It has there- 
fore been asserted that the quantum stress field is not a well-defined physical quantity, even 
though physical intuition may suggest otherwise. A traditional way to develop a quantum 
stress field formalism is to consider the stress field's relationship with the force field. From 
this perspective, the stress field can be defined as any rank-two tensor field whose divergence 
is the force field of the system: 

F a = Vpa ap . (1) 

(Note that the Einstein summation convention for repeated indices is used throughout the 
Letter.) One can add to a af3 a gauge of the form 

^7^ 7 (x), (2) 

where A al31 is any tensor field antisymmetric in (3 and 7, and recover the same force field, 
thereby demonstrating the non-uniqueness of this stress field definition. General formu- 
lations for computing non-gauge-invariant stress fields in quantum many-body systems 



have been derived by Nielsen and Martin, Folland, Ziesche and co-workers, and Godfrey 
P^,P!7| -P^[. There have been several attempts to overcome this problem of non- uniqueness. 
For example, the stress field formalism of Chen and co-workers has been applied to numer- 
ous solid state systems to determine the local pressure around a region . However, their 
method assumes that the potential is pair-wise only. Several ab-initio quantum stress field 
formulations have been developed, as well. Ramer and co-workers developed a method to 
calculate the resultant stress field from an induced homogeneous strain They incorpo- 
rate the additional constraint that the field must be the smoothest fit to the ionic forces. 
This method cannot be used to calculate the residual stress field at equilibrium, nor can it 
determine the energy dependence on strains which do not have the periodicity of the unit 
cell. Filippetti and Fiorentini developed a formulation of the stress field based on the energy 
density formalism of Chetty and Martin p2|j23|| . Since this formulation is based explicitly on 



the energy density, which is not gauge-invariant, the resultant stress field is not unique. Mis- 
tura succeeded in developing a general gauge-invariant formalism for pressure tensor fields 
of inhomogeneous fluids within classical statistical models using a Riemannian geometric 
approach p3|] . 

This Letter extends Mistura's work, developing a Riemannian geometric formalism for 
computing gauge-invariant stress fields in quantum systems within the local density approx- 
imation (LDA) of density functional theory (DFT). We show that the response of the total 
ground-state energy of a quantum system to a local spatially varying strain is a unique and 
physically meaningful field quantity which can be determined at every point in the system. 

Using a procedure well known in continuum theory (see for example Ref. [p5|), one 



can formally relate a Riemannian metric tensor field g a p{^) with the strain field e a/ g(x): 
9ap = <W + 2e Q/3 , with 

e«/3 = ^ {ba-fdpu* + <5 7/3 <9 a u 7 + 5 1K d ot u 1 dpu K ), (3) 



where d a = d/dx a |26|| . Here the strain field is defined in terms of a vector displacement field 
u a which maps coordinates x a in the non-deformed system, to the coordinates x' a = x a + u a 



in the deformed system. 

The stress field a a ^ (x) and strain field obey a virtual work theorem expressing the energy 
response to variations in the strain: 

6E = J ^ga al3 5e aP d 3 x, (4) 

where g(x) is the determinant of g a p(x). It can be shown that the stress field is related to 



the functional derivative of the energy with respect to the metric field |p4 

1 8E 2 5E 



a af) 



^8e aP y/g 5g a p ' 
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We now derive the quantum stress field of a many-electron system in the presence of a 
fixed set of classical positive charged ions using local density functional theory [p7jp8|. The 
ground state electronic charge density of the system is written as n(x) = J2i 4>l ( X )<M X ); 
where <pi are single-particle orthonormal wavef unctions. For this derivation, we assume 
orbitals with fixed integer occupation numbers. The extension to metals with Fermi fillings 
is straightforward, simply necessitating use of the Mermin functional instead of the total 



energy [29|. The total charge density of the system can be written as a sum over all ionic 



charges and n: 



P(*) = E^(x-R*)-^(x), (6) 

i yd 



where Zi is the charge of the i-th ion located at position Rj, and the presence of ^fg insures 
proper normalization of the delta function. The energy of the system can be written as the 
following constrained functional: 

E = E k + ^coulomb + £ xc - E X i (/ V9<P*<Pid 3 x - l) . (7) 

Here E^ is the single particle kinetic energy, ^coulomb is the classical Coulomb interaction 
between the total charge density and itself, and E xc is the exchange-correlation energy of the 
electrons. The appearance of the last term in Eq. [7| is due to the orthonormality constraint 



of the orbitals. (We choose a unitary transformation on {(pi} which enforces orthogonality.) 
One can express E as an integral over an energy density ||23|| . The choice of energy density 
gauge will not affect the derived stress field, since all our results depend only on the total 
energy E. For convenience, we express the energy terms in Eq. [7| as the following: 

^Coulomb = J y/g (pV - ^^^afp) d 3 X 

Exc = J \fg ne hDA (n)d 3 x, (8) 

where T a = —d a V is the electric field due to the Coulomb potential V generated by p, and 
^lda(^) is the LDA exchange-correlation energy density. To obtain the electronic ground- 
state energy, we require 8E /8(f)* = with the additional constraints of a fixed metric (5g a p = 
0) and a fixed ionic charge density (8p = —5n). This implies that the orbitals must obey 
the Euler-Lagrange equations 

- ±f a (Vg~g a %^) + 5 -^^(p l + = Xi(p h (9) 

which can be considered the Kohn-Sham equations in curvilinear coordinates. Also, a least- 
action principle for -E'couiomb requires that p and V obey the Poisson equation: 

^pd a (Vgg a %v) = -An p. (10) 

We now vary the total energy with respect to the metric. It can be proven that we do not 
need to consider variations in the electronic wavef unctions, charge density and potentials, 
since all such variations would vanish due to Eq. |9| and Eq. [10]. This is the same principle 
used in the derivation of the Hellmann-Feynman force theorem and the energy-momentum 
tensor (the variation of the action with respect to metric) in general relativity [pO|-|32|. If a 
different electromagnetic gauge is chosen, variations with respect to the potential V would 
have to be computed explicitly; this change has no impact on the stress field in Eq. [TT| . 
Performing the variation of the total energy with respect to the metric gives the stress field 
in local density functional theory as: 
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v«P = - E daffidpfa + ^-T a Ta + g a p ( - Yl fyfiVfa 

i \ i 

~ g^r^ + - £ 0*0, [Aj + V^j , (11) 

where we have used the relation dy/g/dg alS = —\-Jgg a p- Using Eq. |9|, we can rewrite Eq. 
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as 



(Tap = ~ Yl 9 a(t>*d/3<f>i + J-^a^p + 9ap ( ~ J2 dyftfPfc 
i \ i 

+ 2^ E iVgg^i) - ^T.r + n ^ LD A(n) - S -^j y (12) 

When Eq. [L2] is evaluated at the Euclidean metric (g a p = 5 a p), it gives the stress field at 
zero applied strain, and the {0,} are then solutions to the standard Kohn-Sham equations. 
From here on, we will refer to o a p with an implied evaluation at the Euclidean metric. 

It is important to note several key features of the form of a a p. First, the Coulombic 
contribution to the quantum stress field is equivalent to the classical Maxwell stress field. 



This Coulombic term can be obtained by Filippetti and Fiorentini's formalism in Ref. |22 
if one chooses the Maxwell gauge. Also, the contribution of the exchange- correlation energy 
to our stress field is only in the diagonal (pressure-like) terms, which is the proper behavior 



for local density functionals [14] and is identical to the exchange-correlation stress derived 



in Ref. ||22|| . However, the kinetic contribution to Eq. [12] is unique to our derivation. Our 



stress field contains diagonal terms which are similar to the symmetric and antisymmetric 
kinetic energy densities. By integrating the stress field over all space, we can obtain the 
total stress T a p: 

T a p= 1 1 - V d a (j)*dp4>i + —T a Tp - Sap—F^J 71 

+ 5 a/3 n ^ldaW - -J^j | d 3 x, (13) 

which is identical to the expression derived by Nielsen and Martin [[TJ[ . 

In order to demonstrate the utility of our stress field formalism, we compute within DFT 
the pressure field, | (an + 022 + 033), for two phases of solid molecular hydrogen under exter- 



nal hydrostatic pressure of 50 GPa. [Q. Both structures consist of stacked two-dimensional 



triangular lattices of hydrogen molecules, with the molecular axis parallel to the stacking 
direction and a repeat unit of two layers. The m-hcp structure has alternating layers shifted 
so that each hydrogen molecule is directly above triangular hollow sites in the neighboring 
layers. The second structure (belonging to the Cmca space group) has a different shift, so 
that each molecule lies directly above midpoints between nearest-neighbor pairs of molecules 
in adjacent layers |34|| . The energetics and electronic properties of both structures have been 



studied extensively from first principles [pq , pq] . 

Examination of the pressure field permits us to rationalize the energy ordering of these 
structures. The m-hcp structure is energetically favored by 60 meV/molecule. Figure |l]A 
shows a contour plot of the pressure field for the Cmca structure. The pressure is tensile 
(greater than zero) through the interstitial region, indicating that contraction is locally fa- 
vorable. The pressure is greatest in the volume directly above and below each molecule, 
averaging 3 eV/A 3 . This implies that the system would energetically favor increased inter- 
molecular coordination. In Figure |1]B we show a similar plot for the m-hcp structure. Again 
the pressure within the interstitial region is tensile. However, the pressure field has signif- 
icantly rearranged, and the pressures in the regions above and below molecules have been 
reduced to approximately 2.25 eV/A 3 . It is also clear from the charge density plots (Figure 
|l]C and |l]D) that the reduction in pressure is correlated with an increase in bonding between 
molecules and increased charge derealization. The pressure fields within the molecules are 
compressive, and they are several orders of magnitude larger than the interstitial features. 
However, because these large fields are very similar in both structural phases (and the free 
molecule), they are not important for understanding relative phase stability. Thus, changes 
in the pressure stress field highlight regions and charge density features that contribute to 
favorable energetic changes. 

We have developed a formulation for unambiguously determining the stress field in an 
interacting quantum system described by local density functional theory. The resultant 
expression for the stress field is gauge-invariant with respect to choice of energy density and 
can be obtained via a method analogous to the computation of Hellmann-Feynman forces. 
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Our application of this formalism to solid molecular hydrogen demonstrates that a stress 
field analysis can associate energetics with particular microscopic structural features of a 
material. This stress field formulation has the potential to become an invaluable aid to the 
understanding of structural phenomena in complex solid-state systems. 
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FIGURES 



FIG. 1. Contour plots of the pressure field and charge density within DFT for the Cmca 
structure (panels A and C) and for the m-hcp structure (B and D) of solid hydrogen. The vertical 
axis is the stacking direction for the layers, and the horizontal axis is the direction along which 
alternating layers are shifted. The plots are 6.794 A high and 7.206 A wide. Ten contours are 
shown, over a range of 0-3 eV/A 3 for the pressure fields, and from 0-2.5 e/A 3 for the charge 
densities. 
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